Noise, Synchrony and Correlations at the Edge of Chaos 
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We study the effect of a weak random additive noise in a linear chain of N locally-coupled logistic 
maps at the edge of chaos. Maps tend to synchronize for a strong enough coupling, but if a 
weak noise is added, very intermittent fluctuations in the returns time series are observed. This 
intermittency tends to disappear when noise is increased. Considering the pdfs of the returns, we 
observe the emergence of fat tails which can be satisfactorily reproduced by q-Gaussians curves 
typical of nonextensive statistical mechanics. Interoccurrence times of these extreme events are also 
studied in detail. Similarities with recent analysis of financial data are also discussed. 

PACS numbers: 74.40.De, 05.45.Ra, 87.19.1m 



Since the nonlinear phenomenon of synchronization 
was first observed and discussed in the 17th century by 
Huygens, it has become of fundamental importance in 
various fields of science and engineering. It is frequently 
observed in complex systems such as biological ones, or 
single cells, physiological systems, organisms and even 
populations. Synchrony among coupled units has been 
extensively studied in the past decades providing impor- 
tant insights on the mechanisms that generate such an 
emergent collective behavior (H-Q- In this context cou- 
pled maps have often been used as a theoretical model @ . 
Actually, many biological complex systems operate in a 
noisy environment and most likely at the edge of chaos 
@, H[ ■ Therefore studying the effect of a weak noise in 
this kind of coupled systems could be relevant in order to 
understand the way in which interacting units behave in 
real complex systems like for example living cells d, [l(| ■ 
Generally speaking, random noise is considered a distur- 
bance, i.e. something to avoid in order to obtain precise 
measurements or to minimize numerical errors. But, if 
on one hand it is very difficult to completely eliminate 
the effect of noise, on the other hand it can frequently 
have even a beneficial role. Among the many examples 
in physics and biology we may cite stochastic resonance 
[Tlj , noise enhanced stability [HI, induced second-order 
like phase transitions 
munication networks 




or enhanced diffusion in corn- 
Recently random strategies 
have been demonstrated to be very successful also in mi- 
nority and Parrondo games [lol [lrjj and in sociophysics 
models related to efficiency in hierarchical organizations 
[l7l [HI or even in Parliament models [l9| . 



In previous studies [20, [21| the effect of a small noise 
on globally coupled chaotic units was presented for sev- 
eral kind of systems and a universal behavior related to 
the Lyapunov spectrum was found to be a common fea- 
ture. Power-law correlations and intermittent behavior 
have also been observed in lattices of logistic maps when 



FIG. 1: (Color online) A pictorial view of a chain of N = 50 lo- 
cally coupled logistic maps embedded in a noisy environment. 
The different colors indicate, at a fixed time t, different values 
of the maps in the interval [—1,1]. 



some kind of global coupling exist among them, see for 
example refs.JH, H3- I n this letter, we consider a linear 
chain of N locally coupled logistic maps and we explore 
the role that a small random noise can have in creating a 
strong intermittent behavior and its influence on the syn- 
chronization patterns. We consider only local coupling 
and long-range correlations induced by the noisy envi- 
ronment in which our maps are embedded. At variance 
with previous studies, maps are not in a chaotic regime, 
but at the edge of chaos, where the Lyapunov exponent 
is vanishing (2^ [HI- Moreover, in order to investigate 
our intermittent behavior, we study the pdfs of the re- 
turns of our fluctuating time series, as successfully done 
in different contexts for models showing Self-Organized 
Criticality [26l - [28| . In this respect our model presents 
new features never addressed before as far as we know 
and more related to realistic complex systems. Our re- 
sults can also be framed in the context of non extensive 
statistical mechanics [29U31I ] and analogies with recent 
findings for stock market data analysis [HI, [33[ will be 
addressed. 
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FIG. 2: We show the effect of noise in the normalized returns 
of Eq. (2) for the case N = 100, fi = fi c = 1.4011551..., 
e = 0.8 and r = 32 time steps. At time t = 15000 we switch 
on the noise, with a m ax = 0.002 in panel (a) and o max = 0.2 
in panel (b), then we follow the maps for 100000 iterations. 
See text for further details. 



The model of a linear chain of N coupled logistic maps is 
the following 



= (1 - e) / (xt*) +-[f (xt* 



l )+/0* m )]+°-(t) 

where e e [0, 1] is the strength of the local coupling of 
each map with its first neighbors sites on the chain and 
the additive noise a(t) is a random variable, fluctuating 
in time but equal for all the maps, uniformly extracted in 
the range [0, u max \- In our case the i-th logistic map at 
time t is in the form / (xf) = 1 — fx (xj) , with /i £ [0, 2] 
and with / (xlj taken in module 1 with sign (in order 
to fold the maps' outputs back into the [—1,1] interval 
when the noise takes them out of it). The system has 
periodic boundary conditions. See Fig.l for a pictorial 
view. 

In the absence of noise, this model was extensively stud- 
ied by Kaneko et al. [6|, in particular in the chaotic 
regime, where the coupled maps show different patterns 
of synchronization as function of the coupling strength 
e. Here we consider the effect of a variable addition of 
noise on the same system, but at the edge of chaos, i.e. 
at fi c = 1.4011551... . Following a procedure adopted 
in ref.[23], in order to subtract the synchronized compo- 
nent and keep the desynchronized part of each map we 
consider, at every time step, the difference between the 
average < x\ > and the single map value x\. 



Then we further consider the average of the absolute 
values of these differences over the whole system, i.e. 
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FIG. 3: Asymptotic pdfs of the returns for N = 100 maps 
at the edge of chaos, i.e. /i c = 1.4011551..., with e = 0.8, 
r = 32 (after 100000 iterations) and for increasing values of 
the noise: cr max = 0.002 (a), a max = 0.01 (b), a ma x = 0.05 
(c) and <Jmax = 0.3 (d). Fat tails are more evident for weak 
noise and tend to diminish by increasing noise. We report 
also fits of the simulation data (full curve) by means of q- 
Gaussian curves with values q = 1.54, q = 1.41, q = 1.1 and 
q = 1 (corresponding to a Gaussian) respectively. Returns are 
also normalized to the standard deviation in order to have pdf 
with unit variance. See text for further details. 



d t = j^T,^L 1 \xl— < x\ > |, in order to measure the 
distance from the synchronization regime at time t. If 
all maps are trapped in some synchronized pattern then 
this quantity remains close to zero, otherwise oscillations 
are found. As commonly used in turbulence or in finance 
[2fjl |32| . , we analyze these oscillations by considering 
the two-time returns Adt with an interval of r time 
steps, defined as Ad t = d t+T — d t . 

In Fig. 2 we show that this quantity is very sensitive 
to the noise intensity. In panels (a) and (b) we plot 
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FIG. 4: Top panel: The q-logarithm of the pdf reported in 
Fig. 3 (a) (normalized to the peak) is plotted as function of x 2 . 
A q- logarithmic curve with q — 1.54 fits the points with a 
correlation coefficient equal to 0.9958. Bottom panel: Same 
simulation reported in Fig. 3(a), but with the maps in the fully 
chaotic regime (/i = 2). See text for further details. 



the time evolution of Adt (normalized to the standard 
deviation of the overall sequence) for two different 
simulations obtained with a linear chain of N = 100 
maps, with e = 0.8, r = 32 and considering the maps at 
the edge of chaos. For both the simulations we consider 
a transient of 15.000 iterations, during which the system 
evolves in the absence of noise (cr max = 0), then we 
suddenly increase the level of noise bringing it on at 
Pmax = 0.002 (a) and a max = 0.2 (b) respectively: 
it clearly appears that only in presence of weak noise 
(panel (a)) the returns time series shows large deviations 
from the synchronized pattern of the transient, while 
a higher noise intensity destroys the intermittency and 
induces Gaussian fluctuations. 



In order to quantify such a different noise-dependent be- 
havior we plot in Fig. 3 the asymptotic probability distri- 
bution function (pdf) of normalized returns for increas- 
ing values of a max , from 0.002 to 0.3. Fat tails in the 
pdfs are visible only for a max < 0.1 and, following what 
was already done for a single logistic map at the edge of 
chaos (23.125}. we tried to reproduce them by q-Gaussian 
curves, typical of non extensive statistical mechanics [29[ 
and usually found in complex systems presenting var- 
ious kinds of correlations. q-Gaussians are defined as 

1 

G q (x) = A [l — (1 — q) (3x 2 ] 1 ~ q , where q is the entropic 
index (which evaluates deviations from Gaussian behav- 
ior), 1//3 plays the role of a variance and A is a normaliza- 
tion parameter. These curves actually fit very well the 
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FIG. 5: Asymptotic pdfs of the returns for N = 100 maps 
at the edge of chaos, i.e. /i c = 1.4011551..., with e = 0.8, 
t — 32 as in the case of Fig. 3(a), but now taken with higher 
statistics in order to test the robustness of previous results. 
In this case we considered 10 6 iterations and two histograms 
with a different number of bins, i.e. (a) 150 and (b) 400 re- 
spectively. The same q-Gaussian curve reported in Fig. 3(a), 
here shown for comparison, reproduces very well the new nu- 
merical simulations independently of the size of the bins used 
for the histogram. 



numerical pdfs, as also reported in the panels of Fig. 3 
(full lines). For a max — 0.002 (panel (a)), where the tails 
are very pronounced, one has q ~ 1.5 while, for higher 
values of noise, the tails tend to disappear and the value 
of q decreases asymptotically towards q = 1, which cor- 
responds to a Gaussian pdf (panel (d), with a max = 0.3). 
This definitively demonstrates that if some noise creates 
intermittency and correlations, too much noise destroys 
them. 

As further test to verify the accuracy of the q-Gaussian 
fit shown in Fig. 3 (a), in the top panel of Fig. 4 we plot 
(as open circles) the g-logarithm (defined as ln g z = 
\z x ~ q — 1]/[1 — q], with lni z = In z) of the corresponding 
pdf, normalized to its peak, as function of x 2 , and we 
verify that a y-logarithm curve with q = 1.54 fits very 
well the simulation points with a correlation coefficient 
equal to 0.9958. On the other hand, in the bottom panel 
of Fig. 4, we show that the Gaussian behavior of Fig. 3(d) 
can be also obtained considering the same parameters of 
Fig. 3 (a), but with the maps in the fully chaotic regime, 
i.e. with fj, — 2 instead of fi = \i c . This indicates that the 
edge of chaos condition is strictly necessary for the emer- 
gence of intermittency and strong correlations in pres- 
ence of a small level of noise. Of course the cases fi = fi c 
and fi = 2 are two limiting ones. We also checked that, 
changing the order parameter in the interval \p c , 2] , other 
kinds of non-gaussian pdfs occurs, which are very often 
asymmetric or a superposition of these two extreme cases. 
Finally, in order to test the robustness of the fat tails pre- 
viously found, we show in Fig. 5 the same case reported in 
Fig. 3(a) but now obtained with more statistics, i.e. con- 
sidering 10 6 iterations and pdfs calculated with different 
number of bins. The same g-Gaussian curve of Fig. 3(a) 
(full line) continues to reproduce very well the new data 
both in the central part and in the tails. 
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FIG. 6: A summary of the main results found in our study 
as a function, respectively, of: the size of the system N (a), 
the level of noise a max (b) , the coupling e (c) and the returns 
interval r (d). The maps were always considered at the edge of 
chaos (where the system has a multifractal structure which is 
consistent with the r- independence of q) . See text for further 
details. 
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FIG. 7: Left column panels: plots of the interoccurrence times 
Tj of returns for increasing thresholds in the case N = 100, 
e = 0.8, (Tmax = 0.002 and r = 32. Right column panels: Pdf 
of the time series reported on the left panels. These pdfs are 
nicely fitted by gi-exponential curves, whose value of qi is also 
reported. See text for further details. 



If one considers the value of the entropic index q, 
emerging through g-Gaussian fits of the returns Pdfs, as a 
measure of the correlations induced by the noisy environ- 
ment on our chain of coupled maps at the edge of chaos, 
it is worthwhile to explore how this value changes as func- 
tion, not only of the noise a max , but also of the number 
N of maps, the coupling strength e and the returns time 
interval r. We show in Fig. 6 a summary of the results 
obtained in this direction for fj, = fj, c = 1.4011551... and 
changing the parameters a max = 0.002, N = 100, e = 0.8 
and t — 32 one at a time and then calculating the corre- 
sponding values of q as reported. More precisely, in panel 
(a) we plot the entropic index as function of 1 /N and we 
see that q remains greater than 1 also for very large N, 
thus implying that the noise induced correlations are not 
a finite-size effect. The influence of noise on the value of q 
used to fit the pdf of the returns is reported in panel (b) , 
where an asymptotic convergence towards 1, for strong 
noise, and towards ~ 1.5, for weak noise, is clearly visi- 
ble. Quite interestingly, the plot of q as function of the 
coupling strength, panel (c), has a maximum in corre- 
spondence of e ~ 0.8, a value which evidently allows an 
optimal spreading of correlations over the maps chain in 
the presence of a small noise and corresponds to the fat- 
test tails in the Pdfs, i.e. to frequent large jumps. The 
fact that they occur for finite levels of e is somewhat rem- 
iniscent of phenomena such as stochastic resonance fTJlj j . 
Finally, in panel (d) we plot q versus the time interval r 
used to calculate the returns: the resulting points seem 
to stay constant in the range of r explored, thus con- 
firming the robustness of this kind of correlations for low 
levels of noise. 
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FIG. 8: Panel (a): The values of the index qi, resulting by 
(^-exponential fits of the interoccurrence time series Pdfs, are 
reported in correspondence of five values of the threshold (1.5, 
1.75, 2.0, 2.25, 2.5). A linear fit is also plotted for compar- 
ison. Panels (b): The auto-correlation function Ct/,(A) for 
the interoccurrence time series is plotted for the three corre- 
spondent thresholds values of Fig. 7 and the numerical points 
are fitted by power-law curves Cth(A) ~ A~ 7(t,l) with "f(th) 
equal to, respectively, 0.42 (th = 1.5), 0.37 (th = 2.0), 0.34 
(th = 2.5). 



Long-term correlations in a system typically yield power- 
law asymptotic behaviors in various physically relevant 
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properties. In studies of financial markets, it was re- 
cently observed [32| power-law decays in the so-called 
'interoccurrence times' between sub sequential peaks in 
the fluctuating time series of returns like those shown in 
Fig. 2(a). If we fix a given threshold, the sequence of the 
interoccurrence time intervals results to be well defined 
and it is then possible to study its pdf. We do this for our 
usual chain of N = 100 maps at the edge of chaos, with 
e — 8, r = 32 and for a weak noise with a max = 0.002. 
In the left panels of Fig. 7, the interoccurrence time series 
for the normalized returns are plotted (from top to bot- 
tom) in correspondence of three increasing values of the 
threshold (1.5, 2.0 and 2.5), while the correspondent pdfs 
are reported on the right. In all the cases c^-exponentials 
(i.e., Pdf oc [1 — (1 — (J^Tj/VgJ 1 / 1-9 *, where the subindex 
i stands for interoccurrence) satisfactorily fit the data 
for values of qi which depend on the threshold, in com- 
plete analogy with what was observed for financial data 

In Fig. 8(a) we also show that qi scales linearly as func- 
tion of the threshold. This can be considered as a 
further footprint of the complex emergent behavior in- 
duced on the system by the small level of noise con- 
sidered. Interestingly enough, in the limit of van- 
ishing threshold, q^ approaches unity, i.e., the behav- 
ior becomes exponential, which is precisely what was 
systematically observed in financial data [33[. Fi- 
nally, we also calculated the auto-correlation function 

C th {A) = A'J2 L k ~ A (n(k)- < r % >)( n (k + A)- < n >) 

(ACF) for the interoccurrence time series reported in the 
left panels of Fig. 7, where L is the length of the time 
series, th stands for threshold and A 1 is a normalization 



factor. As shown in Fig. 8(b), for the corresponding val- 
ues of the threshold considered, we found a power-law 
decay Ct/,(A) ~ A -7 ^) with values for the exponent 
^(th) decreasing with the increase of the threshold and 
included in the interval [0.34,0.42], in a gree ment with 
analogous results found in financial data [32l |33| • This 
shows also the presence of memory effects induced by 
noise, in addition to the correlations already pointed out 
by the deviations from Gaussian behavior quantified by 
the entropic index q. 

In conclusion, we have studied the effect of a small ad- 
ditive noise on a synchronized linear chain of N locally- 
coupled logistic maps at the edge of chaos. We found 
strong intermittent fluctuations in the returns, whose 
pdfs arc well fitted with g-Gaussians. The corresponding 
interoccurrence times for fixed threshold exhibit strong 
analogies with financial data. This behavior could bring 
interesting insights on the several common features of real 
systems of different nature which often operate at the 
edge of chaos and in weakly noisy environments. The 
study of further details of this phenomenon in various 
complex systems, including earthquakes, is in progress 
and will be reported elsewhere. 
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